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Abstract 

Recombination  of  a  surface-radical  with  a  surface-olefin  (SR-SOR)  to  form  a  six-membered  ring  is  a  critical  step  in  the 
current  mechanism  for  chemical  vapor  deposition  growth  of  the  diamond  (100)  surface.  We  estimate  the  rate  constant  for 
SR-SOR  by  combining  quantum  chemistry  calculations,  molecular  mechanics  calculations,  and  transition  state  theory.  The 
ab  initio  calculations  include  extensive  electron  correlation  (MP2  and  GVB  *  SD  Cl)  on  cluster  models  which  were  corrected 
for  steric  interactions  of  the  cluster  with  the  rest  of  the  surface  and  for  strain  effects  on  the  lattice.  The  ab  initio  vibrational 
frequencies,  which  were  used  to  construct  a  partition  function  for  calculating  the  entropy,  were  calculated  at  the  MP2  level. 
Transition  state  theory  was  used  to  obtain  the  rate  constant,  k  =  5.6  X  1012  e_8800/Rr  s-1.  This  implies  that  under  normal 
growth  conditions  SR-SOR  is  fast  compared  to  competing  gas-surface  reactions. 


1.  Introduction 

Because  of  strong  industrial  interest  in  developing 
low-pressure  technologies  for  synthesizing  diamond, 
there  has  been  considerable  effort  in  elucidating  the 
fundamental  mechanisms  of  diamond  film  growth 
[1-3].  Since  the  kinetics  and  thermochemistry  of 
carbon-hydrogen  systems  are  well  known,  the  chem¬ 
ical  vapor  deposition  (CVD)  of  diamond  is  probably 
the  best  candidate  for  understanding  CVD  at  a  de¬ 
tailed  level. 

A  number  of  experiments  have  shown  that  methyl 
radicals  are  generally  the  dominant  gas  phase  precur- 
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sor  species  reacting  with  the  surface  to  grow  dia¬ 
mond  [4-8],  and  several  detailed  chemical  kinetics 
mechanisms  have  been  proposed  to  analyze  growth 
on  various  idealized  diamond  surfaces  [9-16].  The 
overall  process  is  given  by 

H2  -*  2H  ',  (1) 

CH4  +  H'-»CH3  +  H2,  (2) 
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For  tetrahedral  electronic  materials  (Si,  GaAs, . . . ) 
the  (100)  surface  has  generally  been  the  surface  of 
interest  for  commercial  growth  of  thin  films  by  CVD 
and  MBE  technology.  Hence  the  mechanisms  for 
diamond  growth  on  this  surface  are  of  interest.  Both 
AFM  and  STM  measurements  have  established  that 
the  stable  surface  is  hydrogen  terminated  with  a 
2X1  reconstruction  as  indicated  in  Fig.  1 .  A  close-up 
of  the  five-membered  rings  (denoted  as  C5)  is  shown 
in  Fig.  2a.  This  C5  ring  must  be  opened  and  con¬ 
verted  to  a  six-membered  ring  (C5  to  C6)  during 
growth.  The  Garrison-Brenner  dimer  mechanism  has 
been  proposed  to  explain  this  process,  and  it  has 
been  analyzed  in  terms  of  known  rates  for  analogous 
gas  phase  hydrocarbon  reactions.  One  step  in  this 
reaction,  namely  the  surface-radical-surface-olefin 
recombination  (SR-SOR)  step,  involves  surface  con¬ 
straints  that  preclude  using  estimates  from  gas  phase 
reactions.  In  this  Letter  we  report  extensive  ab  initio 
calculations,  including  electron  correlation,  to  exam¬ 
ine  this  step.  We  conclude  that  the  rate  is  fast  enough 


Fig.  1.  Dimer  paired  structure  of  the  H  stabilized  C(100)  surface; 
unshaded  =  H,  shaded  =  bulk  C,  hatched  =  bulk  terminating  C. 


to  be  consistent  with  the  Garrison-Brenner  mecha¬ 
nism. 

In  Section  2  we  review  the  Garrison-Brenner 
mechanism.  Section  3  discusses  the  details  of  the 
calculation,  culminating  in  a  rate  constant  for  the 


(f)  C6*  site 


(e)  C5d*  site 


(d)  C5M*  site 


Fig.  2.  The  Garrison-Brenner  mechanism  for  dimer  ring  opening  during  CVD  growth  of  C(100). 
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SR-SOR  step.  Section  4  discusses  the  relation  to 
other  rates  involved  in  diamond  growth. 


2.  The  mechanism 

The  Garrison-Brenner  [5]  dimer  mechanism 
sketched  in  Fig.  2  was  based  on  molecular  dynamics 
simulations  (using  the  Brenner  hydrocarbon  potential 
[17])  and  has  been  analyzed  in  detail  by  Harris  and 
Goodwin  [15]. 

(i)  It  commences  with  addition  of  a  CH3  radical 
to  a  surface  radical  site  (Fig.  2b-2c),  forming  the 
structure  denoted  as  C5M. 

(ii)  The  next  step  (Fig.  2c-2d)  is  abstraction  of  a 
hydrogen  atom  from  the  CH3  group  to  make  the 
radical  denoted  as  C5M  * . 

(iii)  Through  the  3-scission  electronic  rearrange¬ 
ment  (Fig.  2d-2e),  this  structure  isomerizes  to  an¬ 
other  radical,  denoted  as  C5*d.  Here  the  dimer  bond 
breaks  and  a  double  bond  forms  between  the  carbon 
atom  (Cp)  of  the  original  structure  and  the  carbon 
atom  (Ca)  that  started  as  an  adsorbed  methyl  (now 
the  a-carbon  of  the  olefin).  Such  unimolecular  reac¬ 
tions  are  generally  fast  compared  with  bimolecular 
reactions. 

(iv)  The  final  SR-SOR  step  in  the  mechanism 
(Fig.  2e-2f)  is  intramolecular  attack  of  the  radical 
carbon  (Cr)  of  C5*d  with  the  Ca  of  the  double  bond. 
This  leads  to  a  six-membered  ring  denoted  as  C6’, 
containing  a  radical  site  (Cp).  It  is  this  step  that  we 
examine  herein. 

Gas  phase  reactions  involving  radical  attack  at  a 
doubly  bonded  carbon,  such  as 


H2Ca=CpH2  +  CrH3 - »  H2Ca-CpH2 ,  (5) 

are  typically  very  fast  and  proceed  with  almost  no 
barrier.  However,  on  the  surface  two  factors  consid¬ 
erably  decrease  the  rate: 

(i)  the  initial  equilibrium  distance  between  the 
radical  carbon,  Cr,  and  the  a-carbon,  Ca,  is  2.85  A, 
almost  twice  the  normal  CC  bond  distance,  and 

(ii)  Cr  is  tightly  constrained  by  the  lattice  from 
moving  toward  Ca,  while  moving  Ca  toward  Cr 
bends  the  bonds  Cp  to  the  lattice  or  makes  the  bonds 
at  Cp  nonplanar,  leading  to  strain. 


Thus,  on  the  growing  diamond  surface  the  SR-SOR 
reaction  step 

C5*d  -  C6*  (6) 

could  have  a  significant  barrier,  which  might  make 
the  constrained  SR-SOR  reaction  (6)  sufficiently  slow 
on  the  diamond  surface  as  to  limit  the  rate  of  dia¬ 
mond  formation. 


3.  Calculational  details 

3.1 .  Ab  initio  quantum  chemistry 

In  order  to  obtain  accurate  potential  energy  sur¬ 
faces  for  reactions  such  as  (5)  or  (6),  it  is  essential  to 
include  all  electron  correlation  effects  that  change 
during  the  reaction.  For  this  purpose  we  use  the 
generalized  valence  bond  (GVB)  method  [18],  in 
which  all  electrons  involved  in  bonds  that  change 
during  the  reaction  are  fully  correlated  (these  are  the 
active  electrons),  and  all  other  pairs  of  electrons  are 
calculated  self-consistently  (the  semi-active  elec¬ 
trons).  For  reactions  (5)  and  (6)  there  are  three  such 
active  electrons.  For  the  reactant  these  correspond  to 
the  electrons  of  the  rr-bond  (C„  =Cp)  and  the  radical 
electron  of  Cr.  In  the  product  they  correspond  to  the 
two  electrons  of  the  new  C-C  cr-bond  (Cr-Ca)  and 
the  radical  electron  on  Cp.  There  are  two  ways  to 
spin  pair  these  three  electrons  (one  corresponding  to 
reactant  and  the  other  to  the  product)  and  the  GVB 
method  calculates  all  orbitals  self-consistently  (these 
three  plus  the  orbitals  for  all  other  electrons)  while 
optimizing  the  spin  coupling  [19].  This  type  of  calcu¬ 
lation  is  also  referred  to  as  complete  active  space- 
self-consistent  field  (CASSCF).  After  calculating  the 
GVB  wavefunction,  we  then  allow  all  single  and 
double  excitations  from  the  three  GVB  configura¬ 
tions  to  all  possible  virtual  orbitals,  GVB  *  SD  Cl. 
This  accurately  approximates  the  results  of  very 
complete  Cl  calculations  [20].  A  second  approach 
(MP2)  starts  with  the  Hartree-Fock  wavefunction 
and  uses  second-order  (Moller-Plesset)  perturbation 
theory  to  include  the  dominant  many-body  effects 
[21]. 

To  describe  the  reaction  path  we  considered  40 
positions  of  the  critical  carbon,  C„,  selected  to  de- 
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Fig.  3.  The  energy  surface  for  the  SR-SOR  reaction  from  GVB-SD  Cl  calculations.  The  reactant  site  is  denoted  as  0.0.  The  saddle  point  for 
the  reaction  is  denoted  as  9.3.  The  product  is  far  to  the  left  and  top  of  the  figure  with  an  energy  of  -  35.1  kcal/mol.  The  contour  spacing  is 
1.25  kcal/mol.  The  coordinates  X  and  Y  are  the  Cartesian  coordinates  for  motion  of  CQ  parallel  (X)  and  perpendicular  (y  )  to  the  surface 
while  remaining  in  the  (110)  plane  of  the  crystal  which  contains  Cr  and  Cp.  Tick  marks  are  0.1  and  0.05  A  in  the  a:  and  y  directions, 
respectively. 


scribe  the  reactant,  saddle  point,  and  product  states. 
We  fixed  Cr,  Cp  and  the  bonds  to  Cr  and  on  Cp 
representing  the  bulk  atoms.  Then  for  each  position 
of  Ca,  we  optimized  at  the  MP2  level  the  position  of 
the  three  hydrogens.  For  each  such  geometry  from 
MP2,  we  calculated  the  GVB-CI  and  GVB  *  SD  Cl 
energies. 

The  above  calculations  were  carried  out  using  the 
C3H7  cluster  obtained  from  Fig.  2e  by  using  H’s  to 
represent  the  subsurface  carbons  (keeping  the  same 
angles).  We  also  performed  ab  initio  calculations  on 
the  C10H,5  cluster  shown  in  Fig.  2e  and  obtained 
from  Fig.  1  by  using  hydrogens  to  represent  the 
subsurface  (shaded)  carbons.  This  allowed  us  to 
estimate  the  effects  of  cluster  size  on  the  reaction 
barriers  by  considering  (i)  the  reactant  geometry,  and 
(ii)  the  transition  state  geometry.  In  these  later  calcu¬ 
lations,  the  positions  of  Cr,  Ca,  and  and  the  three 
hydrogen  atoms  bonded  to  these  carbon  atoms  were 
all  optimized,  giving  a  better  model  for  relaxation  of 
the  surface.  The  remaining  seven  carbon  atoms  and 
twelve  hydrogen  atoms  were  held  fixed  to  model  the 
lattice.  This  optimization  was  carried  out  at  the  MP2 
level  of  M0ller-Plesset  perturbation  theory  [21].  Us¬ 
ing  the  MP2  optimized  geometries,  we  carried  out 
the  GVB-CI  calculations  to  obtain  the  barrier. 

For  all  GVB  calculations  on  the  three-carbon 


cluster  we  used  the  Dunning/Huzinaga  double  £ 
basis  set  [22,23]  plus  diffuse  s  and  p  functions 
( l s  =  0.0474,  £ p  =  0.0365)  plus  one  set  of  d  polar¬ 
ization  functions  (  £  =  0.75).  The  triple  zeta  contrac¬ 
tion  of  the  6s  basis  was  used  for  hydrogens.  The 
6-3 1 G  *  *  basis  set  was  used  for  the  MP2  optimiza¬ 
tions  on  the  three-carbon  cluster  while  6-3 1G*  was 
used  for  the  MP2  optimizations  on  the  ten-carbon 
cluster.  MP2  was  performed  using  GAUSSIAN  92 
[24],  The  HF,  GVB-CI,  and  GVB  *  SD  Cl  calcula¬ 
tions  used  the  GVB  [25]  1  and  MOLECULE  SWE¬ 
DEN  suites  of  programs  [27]. 

3.2.  Potential  surface 

Using  the  energies  from  GVB*SD  Cl  and  MP2 
calculations,  we  determined  the  potential  energy  sur¬ 
face  shown  in  Fig.  3  (the  reaction  path  is  shown  by 
the  dashed  curve).  The  calculated  reaction  barrier 
generally  falls  with  increasing  levels  of  correlation. 
Thus, 

E™  =18.1  kcal/mol, 


'  See  Ref.  [26], 
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E^2  =  10.0  kcal/mol, 

C1  =10.5  kcal/mol,  and 
E^b*sdci  =  9.3  kcal/mol. 

Hartree-Fock  (HF)  leads  to  a  barrier  about  a  factor 
of  two  too  high  and  to  a  poor  location  of  the 
transition  state  geometry.  The  Eact  =  9.3  kcal/mol 
from  the  GVB’SD  Cl  is  expected  to  be  about  1 
kcal/mol  too  high  from  residual  errors  due  to  in¬ 
completeness  of  the  basis  set  and  the  Cl  expansion 
(see  Table  1).  Thus  for  the  SR-SOR  process,  we 
obtain  a  corrected  GVB  *  SD  Cl  activation  barrier  of 

E”act  =  8.3  kcal/mol. 

We  calculated  the  zero  point  energy  contribution 
to  the  vibrational  adiabatic  barrier  by  computing  the 
MP2  vibrational  frequencies  at  the  transition  and 
initial  states.  This  adds  0.13  kcal/mol  to  the  activa¬ 
tion  barrier,  leading  to 

^duster  ~  8.4  kcal/mol.  (7) 

3.3.  Molecular  mechanics  calculations 

The  strain  energy  imposed  on  the  lattice  by  the 
reacting  surface  species  was  calculated  with  molecu¬ 
lar  mechanics  using  the  MSXX  many-body  force 
field  (fit  to  the  phonon  dispersion  curves  and  elastic 
constants  of  diamond  [28]).  To  model  the  SR-SOR 
reaction  we  increased  the  cubic  unit  cell  by  a  factor 
of  4  in  the  z  direction  and  by  factors  of  5  in  the  x 
and  y  directions  (leading  to  800  atoms  in  the  unit 
cell).  We  then  cleaved  the  (001)  surface  [28],  leading 
to  a  slab  which  was  then  hydrogenated.  The  cluster 
atoms  optimized  in  the  quantum  chemistry  calcula¬ 
tions  were  held  fixed  at  the  initial,  final,  and  transi¬ 
tion  state  geometries,  while  the  remaining  atoms  of 
the  slab  were  optimized  [28].  Strains  and  steric  inter¬ 
actions  included  in  the  ab  initio  cluster  calculations 
were  not  included  in  the  molecular  mechanics  energy 
to  avoid  double  counting  of  this  portion  of  the 
energy.  Cell  structure  and  atomic  coordinates  were 
updated  at  each  optimization  cycle.  The  atomic  coor¬ 
dinates  were  optimized  2  using  conjugate  gradient 


2  These  calculations  used  POLYGRAF  from  Molecular  Simula¬ 
tions  Inc.  (Burlington,  MA). 


Table  1 

Contributions  to  the  activation  energy  for  the  SR-SOR  step 


Calculation 

Energy  (kcal /mol) 

quantum  chemical  calculations 

GVB  ’  SD  Cl 

9.30 

corrections  (basis  set  and  correlation) 

-1.0 

differential  zero  point  energy 

0.13 

net  electronic  structure  barrier 

8.43 

force  field  calculations 

strain  from  balance  of  surface 

0.73 

van  der  Waals  interactions 

-0.34 

total  £acl 

8.82 

techniques  until  the  RMS  force  per  degree  of  free¬ 
dom  was  less  than  0.01  (kcal/mol)/A.  The  lattice 
strain  energy  (not  including  van  der  Waals  interac¬ 
tions)  and  the  van  der  Waals  energy  were  calculated 
for  the  transition  state  and  the  reactant  state.  These 
energies  do  not  include  the  interactions  already  in¬ 
cluded  in  the  ab  initio  calculations  on  the  cluster. 
The  steric  interactions  with  the  surrounding  lattice 
decreases  in  the  transition  state  relative  to  the  initial 
state.  The  strain  imposed  on  the  lattice  at  the  transi¬ 
tion  state  relative  to  the  reactant  state  is  calculated  to 
add  0.73  kcal/mol  to  the  barrier  while  the  van  der 
Waals  interactions  subtract  0.34  kcal/mol  from  the 
barrier.  Thus  these  effects  increase  the  net  electronic 
structure  barrier  by  0.39  kcal/mol  to  yield  a  net 
barrier  of 

E “rf“B  =  8.8  kcal/mol.  (8) 


3.4.  Transition  state  theory 

We  determined  the  entropy  change  AS*  between 
reactants  and  the  transition  state  from  vibrational 
mode  analyses  that  included  the  effects  of  the  con¬ 
straints  on  the  vibrational  levels  during  the  reaction. 
First,  the  Hessian  (second  derivative  matrix)  was 
calculated  using  MP2  theory  at  each  geometry  of  the 
C3H4  cluster  [24].  Next,  the  3/V  —  6  vibrational  lev¬ 
els  were  calculated  by  diagonalizing  the  (mass- 
weighted)  Hessian  (see  footnote  2).  Finally,  these 
energy  levels  were  used  to  construct  a  partition 
function  from  which  the  entropy  was  calculated  (see 
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footnote  2).  For  the  transition  state  the  imaginary 
vibrational  frequency  was  ignored. 

The  TST  pre-exponential  factor  is  given  by 

ek  T 

A  =  ~^e\p(ASt/R),  (9) 

n 

where  In  e  =  1,  kB  is  the  Boltzmann  constant,  T  is 
the  temperature,  h  is  the  Planck  constant,  R  is  the 
gas  constant,  and  AS*  is  the  change  in  entropy 
between  Cd  and  the  transition  state.  We  estimated 
the  entropies  for  the  initial,  transition,  and  final 
states  of  reaction  (6),  leading  to  A  S*  =  —  2.6 
cal/ mol  K  and 

A  —  5.6  X  1012  (10) 

at  1000  K  from  Eq.  (9). 

Combining  (8)  with  (10)  leads  to  a  total  rate 
constant  for  SR-SOR  of 

*sr-sor  =  5.6X1012  e~8800/Jfr  s_l .  (11) 


4.  Discussion 

The  unimolecular  reaction  (6)  on  the  diamond 
surface  competes  primarily  with  abstraction  and  ad¬ 
dition  reactions  involving  gas  phase  H  atoms.  To 
examine  the  relative  importance  of  these  competing 
gas  phase  reactions,  we  compare  reaction  (6)  with 
reaction  (12), 

Cd*  +H**s->Cd,  (12) 

in  which  the  radical  site  recombines  with  a  gas  phase 
H  atom.  At  1200  K  the  characteristic  time  scale  for 
(6)  is  l/kSR.S0R  =  7X  10“ 12  s.  In  contrast,  assum¬ 
ing  kn  is  in  the  range  [29]  1013  to  1014  cm3/mol  s 
with  a  H  atom  concentration  [29]  of  10“10-10“6 
atoms/cm3,  leads  to  characteristic  time  scales  of 
10“3  to  10“8  s.  Thus  the  SR-SOR  reaction  is  104- 
109  times  faster  than  other  steps  in  the  growth 
process,  indicating  that  it  does  not  affect  the  rate  of 
diamond  growth  in  the  Garrison- Benner  mechanism. 

Reactions  like  SR-SOR  may  also  occur  at  steps  or 
edges  and  can  play  a  role  for  the  (111)  and  recon¬ 
structed  (110)  surfaces.  At  such  sites  the  surface 
constraints  may  place  the  initial  structure  for  the 
SR-SOR  step  either  closer  or  farther  from  the  transi¬ 


tion  state,  resulting  in  changes  in  strain  that  could 
decrease  or  increase  the  barrier  relative  to  that  re¬ 
ported  here. 

It  has  been  speculated  [30]  that  SR-SOR  might 
also  be  a  mechanism  for  carbon  incorporation  in 
GaAs  and  Si.  Because  of  differences  in  stability  of 
various  sp"  structures  [31]  for  C  relative  to  Si,  Ga, 
and  As,  such  case  could  have  dramatically  modified 
barriers. 
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